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Table 1 

Couplings of the AFP action currently under investigation. 
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Figure 4. T^/A for the Wilson (squares) and can- 
didate AFP (octagons) actions. 



Figure 5. Potential V(r)/T vs. rT for the AFP 
action at I3c(Nt = 2) (crosses) and I3c(Nt = 4) 
(octagons). 



5. SUMMARY 

This project is incomplete. At present we have 
one candidate AFP action which costs a factor of 
9 in computer time compared to the Wilson ac- 
tion. It shows asymptotic scaling in its bare cou- 
pling with a small T^/A ratio and its potential 
(measured with naive operators) shows better ro- 
tational invariance than the Wilson action. These 
are however not scaling tests. For the potential 
the AFP action shows qualitatively better scal- 
ing. The only quantitative scaling test we have 
completed so far shows excellent scaling both for 
the AFP and Wilson actions as well. That shows 
the importance of testing the scaling behaviour of 
several different operators. We need other tests of 
scaling and we need to develop the formalism of 
improved operators to test the potential. We are 
testing another RG transformation and its AFP 
action. Our long term goal is to apply our tech- 



niques to full QCD. 
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If ^0 is close to the fixed point action S* than 
Uo is close to {U*} and at leading order 

S*(V)= S*(Uo) + T(Uo,V) 

= S*{Uo) + iS'{V)- SoiUo j). (7) 

This equation can be used to calculate the value 
of the fixed point action on the course configura- 
tion {V} if S* is known on the fine configuration 
{Uo}- Using the parametrization for S* and re- 
peating the above procedure on many {V} con- 
figurations leads to a set of linear equations for 
the coefficients c* . 

On very fine configurations the FP action is 
well described by its quadratic form and can be 
used to evaluate S* (Uo) in Eqn. 7. This in turn 
will give S*(V) on coarser configurations that, if 
necessary, could be used to calculate S* in even 
coarser configurations, etc. 

The procedure could fail if in the original 
minimization was not close to the FP action. We 
found that by carefully adjusting the couplings of 
^0 one can make the difference S* (Uo) — So(Uo), 
that measures this correction, arbitrary small. 

4. THE FEW PARAMETER FP AC- 
TION 

Using 0(50) parameters we could get a good 
parametrization of S* down to correlation length 
^ = 0(1). However this action is too compli- 
cated for numerical simulation. While one can- 
not expect that a few parameter action will cor- 
rectly describe the FP action over a wide range 
of correlation lengths, it is possible to get a good 
fit for the FP action with only a few parame- 
ters if we restrict ourselves to a smaller corre- 
lation length range. For that purpose we gener- 
ated coarse configurations with the Wilson action 
between /3 = 5.0 and 6.0 and fitted S*(V) with 
two operators in 3 representations as described in 
Sect. 3. 3. The coefficients of the approximate FP 
(AFP) action is given in table 1. It is not pos- 
sible to find a FP action using only operators in 
the fundamental representation. 

The AFP fits the minimization data well for 
configurations generated with Wilson /3 > 5.7. 
Below that the quality of the fit deteriorates in- 
dicating that other types of operators are needed 



to describe the FP action on configurations with 
smaller correlation length. We expect the second 
type of block transformation, which gives a more 
local quadratic fixed point, will have better prop- 
erties at small correlation lengths. 

4.1. Properties of the AFP action 

The AFP describes the FP on the critical sur- 
face. The properties of the AFP action at finite 
l3 small correlation length range cannot be pre- 
dicted. A detailed scaling test is needed to show 
the properties of the AFP action. 

We have written programs to simulate actions 
like the ones shown in table 1. We have a ver- 
sion for serial machines which uses a combination 
of overrelaxation and Metropolis updating, and a 
parallelized version using the hybrid Monte Carlo 
algorithm. The time to update a link in the ac- 
tion of Table 1 is about 9 times the update time 
of our Wilson code. 

We have been performing one scaling test, a 
measurement of the torelon mass on lattices of 
constant aspect ratio. We set the scale for the 
scaling test by the finite temperature phase tran- 
sition in finite physical volumes. While one does 
not expect improved asymptotic scaling for the 
AFP action it is nevertheless interesting to check 
the perturbative scaling properties and the value 
of the A parameter of the AFP action. Figure 

4 shows Tc/ A both for the Wilson and AFP ac- 
tions. The AFP action not only shows asymptotic 
scaling from Nt = 2 but the value of its A pa- 
rameter is about a factor of ten larger than the 
Wilson one and therefore much closer in value to 
the continuum A parameters. 

As for the quantity G = ^a(L)/Tc, the AFP 
action appears to scale over the range a > 
l/(2Tc), as shown in Fig. 2 (squares). 

Finally, we have measured the potential from 
the correlation function of Polyakov loops. In Fig. 

5 we show V(r)lTc. The notation is the same as 
for the Wilson action in Fig. 3. Compared to the 
Wilson action the AFP action shows much better 
rotational invariance and, what is more impor- 
tant, it shows better scaling between Nt = 2 and 
4. As we mentioned earlier to see the restoration 
of rotational invariance one needs to use perfect 
operators in the potential measurement. 
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ferent scale 2 block transformation with tunable 
free parameters. Both transformations can be de- 
fined as 

^-I3'S'(V) ^ j ^jj^-i3(S(U) + ^T(U,V)) 

where U is the original, V is the blocked link vari- 
able and T(U, V) is the blocking kernel that de- 
fines the transformation. 

The first transformation we used is a modified 
version of the Swendsen transformation where the 
central link is given a different weight than the 
staples. The other transformation uses an APE 
type smearing to create fuzzy links. The blocked 
link is chosen around the product of two fuzzy 
links with a Boltzman factor as indicated in Eqn. 
2. 

On the critical surface, in the limit /3 ^ cxd 
Eqn. 2 reduces to a saddle point problem giving 

S'(V) = min^u} (S(U) - kT(U, V)) . (3) 

The FP of the transformation is determined by 
the equation 

S*(V) = min^u}(S*(U) - kT(U, V)). (4) 

For smooth configurations this equation can 
be solved perturbatively. At smaller correlation 
lengths the configurations are not smooth any 
more and the solution of the equation has to be 
obtained by numerical minimization. The tun- 
able parameters can be used to find the most local 
fixed point within its class. 

3.3. Parametrization of the fixed point ac- 
tion 

We have used a parametrization of the action 
based on powers of the traces of the loop products 
U(C) = Uc Ufi(fi) where C is an arbitrary closed 
path 

S(U) = \J2 (ci(C)(3-ReTriU(C)) 

c 

+ C2(C)(3- ReTriU(C)f + ... 
+ d2{C){ImTr{U{C)f + •••) • (5) 

Since we are searching for a short range fixed 
point first we restricted the sum in Eqn. 5 to the 
28 loops that are length 10 or less and fit into a 



2^* hypercube. For practical purposes one needs 
to reduce the parameter space even further. For 
numerical simulations we consider the projection 
of the action to a subspace of two operators in 
3 representations. We chose the plaquette and 
the twisted 6 link operator {x,y,z,—x,—y,—z) 
for this restricted parametrization. 

3.4. The fixed point at the quadratic and 
cubic level 

For smooth configurations the RG equation can 
be solved analytically. At the quadratic level 
Eqn. 4 translates into a relation between 2-point 
correlation functions. By repeating the transfor- 
mation the equation can be iterated to find the 
fixed point that is expressed in terms of link vari- 
ables. In terms of the loop parametrization this 
will restrict the coefficients of the linear terms 
ci(C). As there are many more loops than links 
the translation of the quadratic fixed point into 
the loop parametrization is not unique. For ex- 
ample only 7 of the 28 loops in the 2^* hypercube 
are independent at the quadratic level. 

It is possible to go one step further and solve 
the fixed point equation at the cubic level. The 
cubic level does not restrict any of the higher rep- 
resentation terms in the action but constrains the 
ci coefficients of 15 out of the 28 loops. 

Since these calculations are analytic we can eas- 
ily investigate the properties of different renor- 
malization group transformations and find the 
one that gives the most local fixed point action. 
The best FP action is quite local especially for the 
second type of transformation and is dominated 
by the plaquette term in both cases. 

3.5. Numerical minimization 

To go beyond leading order we solved the FP 
equation numerically. First we generate an SU(3) 
configuration {V}. This is our coarse configura- 
tion. Next we choose an action that is suffi- 
ciently simple for minimization but close to the 
fixed point action. The minimization of Eqn. 3 
gives the fine configuration {Uq} that blocks into 
{V} 

S'{V)= minu{So{U)+T{U,V)) (6) 
So{Uo) + T{Uo,V). 
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Figure 2. Scaling test for the Wilson action (oc- 
tagons) and AFP action (squares); Tg is defined 
in constant physical volumes. 



Figure 3. Potential T/(r)/T vs. rT for the Wilson 
action at I3c(Nt = 2) (crosses) and I3c(Nt = 4) 
(octagons). 



that lack of rotational invariance of the direct 
Polyakov loop operators does not mean lack of 
scaling and, alternatively, perfect rotational in- 
variance does not necessarily mean perfect scal- 
ing either. Scaling violation is evident from the 
overall difference of the two potentials. 

3. THE CLASSICAL FIXED POINT AC- 
TION 

3.1. Perfect actions 

Consider an SU(3) gauge theory defined on a 
hypercubic lattice. The partition function is 

Z = J DUe-P'(''\ (1) 

where j3S{U) is some lattice representation of the 
continuum action. It is a combination of the 
products of the U^{n) = e-«f^M(") SU(3) matri- 
ces along arbitrary closed loops in arbitrary rep- 
resentation. 

Under repeated real space renormalization 
group transformations the action moves in a mul- 
tidimensional parameter space. On the critical 
surface it runs into a fixed point (FP). The FP 
of the SU(3) gauge theory is at /3 = cxd and has 



one marginally relevant and infinitely many irrel- 
evant directions. The trajectory which leaves the 
FP along the marginal direction is the renormal- 
ized trajectory, RT. 

The FP and any point on its RT describe per- 
fect actions. There are no lattice artefacts along 
the RT. The FP/RT, while describing perfect ac- 
tions, can be very complicated containing many 
non-negligible interaction terms. For practical 
purposes one needs a sufficiently short ranged 
FP/RT that can be parametrized by a few op- 
erators, at least at correlation lengths that are 
relevant to numerical simulations. 

In the present work we concentrate on locating 
a FP on the critical surface. Since the RT leaves 
the critical surface perpendicular to it the FP de- 
scribes the whole RT close to the critical surface. 
We assume that it stays close at smaller corre- 
lation lengths as well. This is an uncontrolled 
assumption and is potentially the most serious 
source of error for this project. 

3.2. The fixed point action 

The location of the FP is not unique. It de- 
pends on the particular choice of renormalization 
group transformation. We worked with two dif- 
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When we started this project we thought that 
data to show scaling or the lack of it for the Wil- 
son action would be readily available in the liter- 
ature. Unfortunately that is not the case. Scaling 
tests can be very sensitive to finite volume effects 
and the systematic errors associated with infinite 
volume extrapolations make many of the existing 
data unusable. This is especially so for [3 > 6.0. 

We have decided calculate physical observables 
for our scaling test in fixed, finite physical vol- 
umes. This way we can avoid the infinite volume 
extrapolation. For one scaling test we consider 
the quantity G = \J a(L)/L where ulL) is the 
string tension on an volume computed from 
the exponential fall-off of the Polyakov line cor- 
relator (torelon mass). If evaluated in constant 
physical volumes, G is independent of the bare 
coupling in the scaling (continuum) limit. Any 
variation of G is due to lattice artefacts. Constant 
physical volumes can be achieved, for example, by 
evaluating G at the critical coupling of the finite 
temperature phase transition aTc = 1/Nt in vol- 
umes = (r#T)^ where r is some conveniently 
chosen aspect ratio. Figure 1 shows G(L) with 
r = 2.0 as the function of aTc for the Wilson 
action. The plot contains data for Nt between 
2 and 14 spanning the coupling constant range 
from l3 = 5.1 to /3 = 6.4. The torelon masses for 
l3 > 5.9 (Nt > 6) are from published data[3]. The 
large errors are due to the uncertainty of the in- 
finite volume finite temperature phase transition. 
The Wilson action appears to show an approx- 
imately 10% scaling violation between Nt = 2 
and 8. 

To avoid extrapolating /3c to infinite volume 
one can set the scale by the finite temperature 
phase transition at finite, constant physical vol- 
umes NgNT, Ns = t-Nt- Strictly speaking there 
is no phase transition in finite volume and we have 
to define what we consider l3c(Ns). We found 
that the Columbia group's definition [4] of /3c is 
the most reliable on small volumes. It is based on 
monitoring the phase of the Polyakov loop and 
requiring 50-50% occupation time for both the 
deconfined and confined vacua. 

This quantity shows smaller scaling violations 
than the previous one. It is displayed as octagons 
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Figure 1. Scaling test for the Wilson action; Tg is 
defined in infinite spatial volumes. 



in Fig. 2. Surprisingly (to us) the Wilson action 
data for this observable is consistent with scaling 
for Nt > 2. 

Alternatively one could use the potential evalu- 
ated at fixed physical volumes and temperature to 
study scaling. However, the potential is derived 
from a correlator of operators and its behavior is 
sensitive to the properties of operators in a way 
that a mass measurement is not. When scaling 
holds the quark-antiquark potential is rotation- 
ally invariant even at small distances, but that 
might not be evident if one measures an opera- 
tor that creates, instead of one localized fermion, 
some spatially extended quark distribution. The 
naive Polyakov loop is such an operator. The 
perfect potential operator is an operator that de- 
scribes one heavy quark at a definite location. It 
is in principle a combination of Polyakov lines and 
Wilson loops and depends on the action itself. 
Fig. 3 shows V(r)/Tc as the function of rTc for 
the Wilson action at two different values of the 
lattice spacing. Crosses and octagons correspond 
to measurements at (3c{Nt = 2) and (3c{Nt = 4) 
on (2#t)^(1.5#t) lattices respectively. The po- 
tentials, measured with naive Polyakov loops, are 
normalized so V{r = l/Tc) = T^. The lack of 
scaling is evident. We would like to emphasize 
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1. INTRODUCTION 

"Perfect actions" are actions which are defined 
on a lattice but show no lattice artefacts. Any 
action on a renormalized trajectory of any renor- 
malization group transformation is perfect, free 
of lattice artifacts. If one can find such an action 
that, in addition to being perfect, is reasonably 
local as well, than the computational overhead 
associated with simulating the more complex ac- 
tion can be easily compensated by its improved 
scaling behavior. This program was successfully 
carried out for the 2 dimensional 0(3) a model[l]. 

We have been trying to construct and test per- 
fect actions for SU(3) pure gauge theory. Our 
ultimate goal is to find a perfect action for QCD 
with fermions where scaling violations are more 
serious. This talk summarizes some of our results 
to date. 

Lattice regularization introduces various arte- 
facts like scaling violation or rotational sym- 
metry breaking. The distortion caused by lat- 
tice artefacts can depend strongly on the spe- 
cific choice of lattice action. Most numerical cal- 
culations to date use the Wilson plaquette ac- 
tion. This action is the simplest and computa- 
tionally the fastest per site but the continuum 
limit is probably reached only at fairly small lat- 
tice spacing and, consequently, typical pure gauge 
or quenched calculations are done on spatial vol- 
umes 32^ or larger. 
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The alternative to decreasing the lattice spac- 
ing in controlling lattice artefacts is to use an im- 
proved action. Although the history of improved 
actions is long they have had a limited effect on 
actual calculations until now. The reason is that 
the computational overhead of simulating these 
actions was not compensated by the possible im- 
provements obtained in the scaling properties. A 
new action obtained by combining a perturba- 
tively improved action with non-perturbative cor- 
rections has been proposed at this conference [2]. 
The theoretical background of our method is com- 
pletely different from that of [2] . 

Our goal here is to find a systematic way to con- 
struct and test a perfect action which is largely 
free of lattice artifacts. In the following we first 
discuss scaling tests for SU(3) gauge theory. Next 
we describe how a fixed point action can be ob- 
tained and investigate a candidate action. 

2. SCALING TESTS 

Scaling means that all physical dimensional 
quantities show the same functional dependence 
on the gauge coupling. It is different from asymp- 
totic scaling where in addition to scaling we 
require that this functional dependence be de- 
scribed by the 2-loop /3 function. 

Scaling violations can be different for different 
quantities. The ratio of certain observables might 
show scaling at smaller correlation length than 
others. Scaling requires the universal behavior of 
all observables. 



